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ABSTRACT 

This  paper  describes  a  novel  sniper  localization  method 
for  a  network  of  sensors.  This  approach  relies  only  on 
the  time  difference  of  arrival  (TDOA)  between  the 
muzzle  blast  and  shock  wave  from  multiple  single¬ 
sensor  nodes,  relaxing  the  need  for  precise  time 
synchronization  across  the  network.  This  method  is  best 
suited  where  an  array  of  sensors,  on  a  per-node  basis,  is 
not  feasible.  We  provide  results  from  data  collected  in  a 
field. 


1.  Introduction 

There  are  several  commercial  sniper  localization 
systems  by  various  vendors  [1,2].  Vanderbilt  University 
developed  a  soldier  wearable  shooter  localization  system 
[3].  All  these  systems  are  based  on  real-time  operating 
systems  such  as  UNIX,  LINUX,  etc.,  and  have  elaborate 
time  synchronization  mechanism.  Time  synchronization 
allows  all  deployed  sensors  to  share  a  common  time 
reference  so  that  they  can  determine  the  exact  time  of 
arrival  (TO A)  for  both  shock  wave  and  the  muzzle  blast 
for  supersonic  gun  fire.  The  commercial  systems  employ 
array  of  microphones  at  each  location  to  determine  the 
angle  of  arrival  (AOA)  of  the  muzzle  blast  and  the 
shockwave.  With  the  knowledge  of  the  difference 
between  the  muzzle  blast  and  shockwave  arrival  times  it 
is  easy  to  determine  the  sniper  location  and  the  trajectory 
of  the  bullet  [1,3].  However,  the  array  based  sensors  are 
less  suitable  for  man-wearable  systems. 

The  goal  of  this  research  is  to  develop  a  sniper 
localization  system  using  distributed  single  microphone 
sensors  or  PDAs  using  only  time  differences  between 
the  shockwave  and  muzzle  blast  at  each  sensor.  When 
the  sensors  are  distributed,  time  synchronization  among 
the  sensors  is  critical.  If  one  uses  a  non  real  time 
operating  system  such  as  Window  CE,  the 
synchronization  among  the  sensors  is  not  guaranteed.  As 
a  result,  the  TO  A  estimates  of  different  events  could  be 
off  by  1-2  seconds  making  the  localization  impossible. 
Figure  la  shows  the  shockwave  and  muzzle  blast  and 
figure  lb  shows  the  recordings  of  the  same  gun  shot  data 
on  four  different  PDAs  showing  bigger  time  delays  due 
to  lack  of  time  synchronization.  However,  if  the  internal 
clocks  of  the  sensors  are  stable,  one  can  estimate  the 
difference  in  TOA  of  shockwave  and  muzzle  blast.  In 
this  paper,  we  provide  sniper  localization  based  on  the 


time  difference  of  arrival  (TDOA)  of  muzzle  blast  and 
shockwave. 

In  section  2  we  present  the  problem  and  derive  the 
necessary  equations  required  to  determine  the  sniper 
location  using  the  TDOA  of  muzzle  blast  and  shock 
acoustic  wave  at  individual  microphone  sensors.  In 
section  3,  we  present  the  algorithm  and  the  sensor 
localization  results  for  real  data  corresponding  to  a  rifle 
fired  from  a  location.  In  section  4,  we  present  the 
conclusions. 
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Figure  1:  (a)  Shock  wave  and  muzzle  blast  of  a  gun  shot, 
(b)  data  recorded  by  four  PDAs 


2.  Sniper  Localization  Using  TDOA  between  the 
Muzzle  blast  and  Shockwave 

We  assume  that  the  sensors  are  single  microphone 
sensors  capable  of  recording  the  acoustic  signals  due  to 
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rifle  firing.  Individual  microphones  can  detect  the 
acoustic  signals  and  are  capable  of  detecting  the  time  of 
arrival  of  muzzle  blast  and  the  shockwave  and  hence  the 
TDOA  between  the  two.  Even  though  the  sensors  are  not 
synchronized  to  a  single  time  frame,  we  assume  that  the 
TDOA  of  muzzle  blast  and  shockwave  at  each  sensor 
can  be  estimated  accurately.  It  is  well  known  from  the 
ballistic  data  that  the  bullet  looses  its  speed  due  to 
friction  as  it  moves  away  from  the  gun.  In  order  to 
develop  the  theory,  we  first  present  the  sniper 
localization  using  a  constant  velocity  model  for  the 
bullet  and  then  change  the  model  for  the  more  realistic 
case  where  the  velocity  is  not  constant.  In  both  the  cases 
we  assume  that  the  trajectory  of  the  bullet  is  a  straight 
line  -  which  is  a  valid  assumption  for  the  distances  up  to 
300  m  [4]. 

Constant  Velocity  Model:  Figure  2  shows  the  geometry 
of  the  bullet  trajectory  and  the  shockwave  cone.  In 
figure  2,  Z  denotes  the  location  of  the  sniper  and  U  is  the 
unit  vector  in  the  direction  of  the  bullet.  As  the  bullet 
travels  at  super  sonic  speed  the  shockwave  generates  a 
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Figure  2:  Geometry  of  the  bullets  trajectory  and  the 
shockwave  cone 

cone  with  angle  # ,  where  sin  #  =  1  /  m  ,  m  is  the  mach 
number.  The  shockwave  propagates  perpendicular  to  the 
cone  surface  and  reaches  the  sensor  Sk  .  The  point  where 
the  shockwave  radiates  towards  the  sensor  is  denoted 
by  Ak .  By  the  time  the  shockwave  reaches  the  sensor 

Sk ,  the  bullet  has  traveled  from  Ak  to  Ck  and  the  miss 
distance  is  given  by  hk  =||^-^||,  where  ||i?||  denotes 

the  norm  of  the  vector  B  .  Fet  yk  is  the  angle  between 
the  trajectory  of  the  bullet  and  the  line  joining  the  sniper 
location  and  the  sensor  location  Sk  . 

Fet  us  denote  the  time  of  arrival  of  muzzle  blast  and  the 
shockwave  as  Tk  and  tk  respectively, 

ri=M;  =  (i) 


where  v  denotes  the  propagation  velocity  of  the  sound. 
Note  that  in  the  time  the  shockwave  propagates  from  Ak 

to  Sk ,  the  bullet  travels  from  Ak  to  Ck ;  thus  the  bullet 
travels  from  Z  to  Ck  during  the  time  period  tk .  Then 
tk  can  be  re-written  as 
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where  ‘T’  is  the  transpose  and  we  used  the  relationship 
that  || Bk  -Z |  is  nothing  but  the  projection  of  the  vector 

Sk-Z  onto  the  trajectory  of  the  bullet  with  unit  vector 
U .  Using  the  trigonometric  relations  that 
| Bk  -  Z||  =  \\Sk  -  Z||  cos  yk  and  hk  =  || Sk  -  Z||  sin  yk  and 
sin  6  =  1  /  m  we  get 


(sin#  cos  yk  +  cos#  sin  yk} 
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From  equations  (1)  and  (3)  we  find  the  TDOA 
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where  dk=v(Tk-tk )  and.  qk  =[l -s,m(d  +  rk)~\ 

We  use  equations  (5)  and  (4)  to  solve  for  the  sniper 
location  and  the  trajectory  of  the  bullet.  We  now 
consider  the  case  where  the  bullet’s  velocity  changes 
during  the  course  of  its  travel. 


Changing  Velocity  Model:  Figure  3  shows  the 
geometry  of  the  bullet  trajectory.  Notice  that  the  cone 
angles  are  different  at  the  point  Ak  when  the  shockwave 

propagates  to  the  sensor  and  at  the  point  Ck  at  the  time 

the  shockwave  reaches  the  sensor.  This  is  due  to  the 
decreasing  speed  of  the  bullet.  For  the  sake  of  analysis 
we  assume  that  the  bullet  travels  at  average  of  speed  of 

mach  mx  from  Z  to  Bk  and  sin  6X=\  I  mx\  then  the 
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mv 


v 


bullet  travels  from  Bk  to  Ck  at  an  average  speed  of  mach 


m2  and  sin 02  =  l/m2. 


microphone  sensors  distributed  closely  in  a  field  as 
shown  in  figure  4.  The  rifle  is  fired  in  several  directions 
from  two  locations  and  the  data  is  collected  for  process- 
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Figure  3:  Trajectory  of  the  bullet  with  different  cone 
angles 

From  figure  3,  we  find  the  TO  A  of  shockwave  is  given 
by 
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where  p  =  ^s'm2  0^  +  cos2  and  a  =  sin-1  ^s’n  ^ 
Then  the  TDOA  is 

Tk~  h  =1^— 3[[l-/?sin(a  +  y*)] 
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where  wk=\- psm(a  +  yk)  and  dk  =v{Tk  -  tk )  •  Just 
as  in  the  case  of  constant  velocity  model,  the  unit  vector 
in  the  direction  of  bullet  trajectory  is  related  to  7  k  by 


IIs*  -A 
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2.1  Implementation  of  Sniper  Localization  Algorithm 
&  Results 

Data  Collection:  In  order  to  test  the  algorithm,  we 
collected  supersonic  rifle  firing  data  with  eight 


Figure  4:  Data  collection  scenario  with  8  microphones 
and  two  sniper  locations 

ing.  We  used  two  locations  for  the  sniper  which  are 
roughly  60  m  apart  as  shown  in  figure  4.  The  sensors 
were  located  about  250  m  down  range.  A  total  of  8 
sensors  were  used  to  collect  the  acoustic  signatures.  The 
data  is  collected  at  100  KHz  in  order  to  capture  the 
shockwave.  The  acoustic  signatures  collected  were 
processed  to  detect  the  shockwave  and  muzzle  blast.  In 
the  case  of  the  constant  velocity  model,  we  solve 
equations  (4)  and  (5)  or  in  the  case  of  changing  velocity 
model  we  solve  the  equations  (7)  and  (4).  In  order  to 
solve  equation  (7),  we  need  to  compute  the  parameters 
{ v,  w,  a,  yk } ,  V  k  e  {l,  2, •  •  • ,  8} .  The  propagation  velocity 

of  sound  v  is  estimated  using  the  meteorological  data, 
i.e,  the  temperature  using  the  formula 


v  =  331.3J1  + 


m/sec 


273.15 

where  the  temperature  r  in  Celsius.  In  order  to  estimate 
the  parameters  a  and  ft  we  first  need  to  determine  the 

mach  numbers  m1andm2.  From  the  ballistics  of  the 

bullets  used  in  the  rifle,  we  used  the  average  speed  of  the 
bullet  from  the  time  the  bullet  emerges  from  the  muzzle 
and  the  bullets  speed  at  a  distance  of  250  m  for  mx  and 

the  speed  of  the  bullet  at  the  distance  250  m  for  m2 . 
Once  7^  and  are  known,  we  use  (6)  to  determine 
a  and  /? .  That  leaves  us  with  determining  the 
parameters  yk  for  all  k.  However,  this  is  a  difficult  one 
to  estimate  without  the  knowledge  of  the  sniper  location. 
This  is  done  iteratively  with  initially  setting  the  values  of 


yk  =  0.01  for  all  k  and  then  estimating  the  sniper 

location  Z.  In  order  to  estimate  Z  using  equation  (7)  we 
first  linearize  it  resulting  in 
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where  n  is  the  number  of  sensors.  From  equation  (9)  we 
estimate  the  initial  value  of  Z  which  is  then  used  in 
estimating  the  values  of  yk  denoted  by  yk  using 
equation  (7),  that  is, 
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The  next  step  is  to  estimate  the  unit  vector  of  the 
trajectory  U  denoted  by  U  using  linear  equation  (4) 
given  by 
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which  can  be  solved  for  U  using  regression.  We  have 
estimated  the  values  of  Z  ,  and  U  which  can  be  used  to 
re-estimate  the  values  of  yk  using  the  equation  (4),  that 
is, 


yk  =  cos' 


(sk-z)Tu 

k-z| 
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Ideally  if  the  values  of  Z  =  Z  ,  and  U  =  U  the  values  of 
^and  yk  calculated  using  the  equations  (10)  and  (12) 


respectively,  should  be  equal  resulting  in  determining 
the  sniper  location  and  the  bullet’s  trajectory. 

Experimental  Results:  From  the  processed  data,  we 
estimated  the  TOA  of  the  shockwave  and  the  muzzle 
blast  which  gave  us  the  TDOA  between  the  muzzle  blast 
and  shockwave.  We  then  applied  the  algorithm  given 
here  to  estimate  the  location  of  the  sniper.  In  order  to  use 
the  above  algorithm,  we  need  the  information  about  the 
bullet  speed  and  the  propagation  velocity  of  the  sound. 
While  the  propagation  velocity  of  the  sound  can  be 
reasonably  estimated  using  the  meteorological  data  such 
as  temperature,  humidity  and  the  wind  velocity,  the 
bullet  speed  varies  from  bullet  to  bullet  due  to  the 
variations  occurring  in  the  manufacturing  process.  One 
can  estimate  the  bullets  velocity  from  the  N-wave 
generated  by  the  shockwave  [1,  4].  This  requires  the 
knowledge  of  the  diameter  of  the  bullet,  length  of  the 
bullet  and  the  shape  of  the  bullet.  However,  this  does  not 
provide  the  velocity  of  the  bullet  from  the  time  it  was 
fired  from  the  gun  to  the  point  where  the  sensors  are 
located.  In  order  to  overcome  this,  we  allow  the 

algorithm  to  estimate  the  Mach  numbers  mx  and  m2 . 

Better  approach  for  estimation  of  ml  and  m2  is  to  use 

optimization  algorithm  such  as  MATLAB’s 
“fminsearch”.  We  used  fminsearch  to  determine  the 

sniper  location  and  optimize  the  values  of  mx  and 
m2  resulting  in  the  minimization  of  the  sum  of  the 
difference  between  yk  and  yk  for  all  k.  These  algorithms 

are  basically  search  algorithms.  Unfortunately  there  are 
several  local  minima  and  hence  the  sniper  localization 
results  are  sensitive  to  the  initial  estimate  of  the  sniper 
location  used.  In  order  to  overcome  this  difficulty,  the 
fminsearch  is  performed  with  starting  points  in  the  3- 
dimentional  grid  with  X  &  Y  changing  from  -100  to  100 
with  an  increment  of  30  and  Z  changing  from  -10  to  100 
with  increments  of  30.  The  resultant  estimates  of  the 
algorithm  are  averaged  to  determine  the  overall 
estimation  of  the  sniper  location. 

In  Figure  5  the  red  stars  denote  the  location  of  the  single 
microphone  sensors,  the  black  stars  indicate  the 
estimation  of  the  sniper  location  using  the  constant 
velocity  model,  and  the  blue  stars  indicate  the  location 
of  the  sniper  using  the  changing  velocity  model. 


(a)  3D  view  of  estimated  sniper  locations 
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Table  1:  Statistics  of  estimated  sniper  location  by  both 
the  models 


We  are  able  to  localize  the  sniper  position  to  within  60m 
for  the  case  where  the  sniper  fired  the  gun  from  a  berm 
on  the  ground.  These  estimate  appears  to  be  biased 
estimates.  In  the  future  work  we  plan  on  investigating 
the  reasons  for  the  bias. 

3.  Conclusion 

In  this  paper  we  have  presented  a  sniper  localization 
algorithm  using  single  microphone  sensors.  The 
algorithm  is  based  on  the  time  difference  of  arrival  of 
muzzle  blast  and  shock  wave  which  can  be  measured  at 
each  microphone  accurately  without  necessitating  time 
synchronization  among  the  sensors.  This  feature  makes 
the  approach  more  readily  usable  for  man-wearable 
sniper  localization  system. 


(b)  2D  view  of  estimated  sniper  locations  with 
constant  vel.  Model 
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(c)  2D  view  of  estimated  sniper  locations  with 
changing  vel.  model 

Figure  5:  Results  of  the  algorithm 


We  presented  two  models  for  estimation  of  the  sniper 

location.  The  merits  of  the  two  will  be  explored  further 

in  later  work. 
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Table  1  presents  the  average  values  of  X,  Y,  and  Z 
coordinates  and  their  variances  for  the  sniper  location 
using  both  models. 


